The Impact of an Oxygen Dopant in an ideal Bi2Sr2CaCu208+5 Crystal 
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Recent scanning tunneling microscopy studies have shown that local nanoscale pairing inhomogen- 
ities are correlated with interstitial oxygen dopants in Bi2Sr2CaCu2 08+4. Combining electrostatic 
and cluster calculations, in this paper the impact of a dopant on the local Madelung and charge 
transfer energies, magnetic exchange J, Zhang-Rice mobility, and interactions with the lattice is 
investigated. It is found that electrostatic modifications locally increases the charge transfer energy 
and slightly suppresses J. It is further shown that coupling to c-axis phonons is strongly modified 
near the dopant. The combined effects of electrostatic modifications and coupling to the lattice 
yield broadened spectral features, reduced charge gap energies, and a sizable local increase of J. 
This implies a strong local interplay between antiferromagnetism, polarons, and superconducting 
pairing. 



The precise role of the atoms lying off the Cu02 planes 
has been an intriguing puzzle in the cupratcs. While 
largely thought to provide a charge reservoir to dope 
holes into the Cu02 plane, it has become clear that an 
understanding of the pairing mechanism will require ad- 
dressing the large variations in Tc arising from the local 
environment surrounding the Cu02 planes [l[. Empiri- 
cally, the role of the apical or axial orbitals has been a 
vehicle linking Tc to either an effective electron hopping 
t' along diagonal Cu-Cu bondsQ, or to the stability of 
the Zhang-Rice sing let (ZRS)@1. However to date these 
arguments have pointed out possible links but offer little 
microscopic reason for the impact on Tc itself. 

Scanning tunneling microscopy (STM) in 
Bi2Sr2CaCu208+5 (Bi2212) has revealed that nanoscopic 
inhomogeneity is correlated with the location of inter- 
stitial oxygen dopant atoms [1] or the superlattice 
modulation[5|. The location of dopants has been cor- 
related with suppressed peak features with larger gap 
energies in the observed local density of states (LDOS), 
and have been associated with local modification of 
superconducting pairing Q. This suggests a non-trivial 
link between the dopant atoms and the electronic 
properties of the material on a local level. 

In metallic systems, such defects can be effectively 
screened and have little impact on the electrostatics of 
the material. This is in contrast to the cuprates, which 
have poor screening along the c-axis and are unable to 
effectively screen the dopant's impurity charge. As a re- 
sult, these dopants and accompanying strutural changes 
may have a substantial impact on the electrostatic prop- 
erties of the material. This may be reflected in quantities 
such as the charge transfer energy A, effective hoppings 
t, t' , magnetic exchange interaction J, or electron-phonon 
(el-ph) coupling strength A. 

From i — J studies 0, it has been argued that dopants 
give an enhancement of J and thus larger gap energies, 



yet the overall shape of the LDOS suggests that inco- 
herence, giving broad spectral features, is an important 
ingredient to understand local pairing modifications. 

In order to quantitatively address these issues elec- 
trostatic Ewald calculations are performed for Bi2212 
supercells to determine the spatial dependence of the 
Madelung energies around atomic sites in the crystal. It 
is found that while Madelung energies on O and Cu are 
spatially varied on the scale of eVs, these changes largely 
cancel and A is slightly increased near the dopant, yield- 
ing an overall suppression locally to J. This information 
is then combined into exact diagonalization (ED) clus- 
ter studies, yielding effective parameters t, t' , J and J'. 
Large 0(1) changes are found both in t' as well as c— axis 
el-ph coupling A, quantities which are known to strongly 
modify a d^.2_y2 pairing interaction. Finally ED cluster 
studies including c-axis phonons are shown to produce 
broadened spectra, a reduced charge gap, in agreement 
with experiments [4|. As a consequence, a sizable local 
increase of J results due to reduced gap via a gain in 
lattice energy, implying a strong interplay between el-ph 
coupling, and local superconducting pairing and antifer- 
romagnetic correlations. 

Electrostatic calculations using Ewald's method are 
performed on a supercell consisting of 3 x 3 x 1 Bi- 
2212 unit cells and containing 270 atoms. We have ex- 
amined supercells up to four times this size and found 
that the results are not qualitatively different. Each 
unit cell contains 2-primitive cells stacked along the c- 
axis for a total of 4 Cu02 planes, as shown in Fig. 
[T^. Using formal valences for the atoms along with 
the known structural dataQ, the Madelung energies <!> 
obtained are ^apex,piane = 18.48, 10.16 eV for the api- 
cal, planar oxygen sites respectively, and ^cu = —38.22 
eV for the copper site, consistent with values reported 
previously The charge transfer A is related to the 
difference in Madelung energies for the Cu and O sites. 
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FIG. 1: (a) Schematic location of the dopant oxygen in Bi-2212 unit cell. Individual Cu02 planes are labeled by the symbols 
shown. The arrows on the left indicate the orientation of the local oxygen crystal fields for the undoped lattice, (b - d) Madelung 
energies of the apical O, planar O and Cu sites, respectively, in the doped lattice. The local value of the charge transfer energy 
A (e) and derived exchange interaction J (f), obtained by averaging A<I>a/ between Cu and its four neighbouring O sites, (g) 
Resulting J for different distributions of charges. All values have been normalized to the undoped lattice values. The distance 
to dopant is defined as the distance to the closest dopant, accounting for the periodicity of the superlattice. 



A$M = $0 — ^Cu and is given by: 
A<I>i\f 

A^^-Icu{2) + Ao{2)~- (1) 
e(oo) dp 

where Icu{2) and Ao{2) are the second ionization and 
electron affinity energies for the Cu and O sites re- 
spectively. The factor of e^/d represents the contribu- 
tion of the Coulomb interaction between the introduced 
electron- hole pair. In this work we take the dielectric 
constant e{oo) = 3.5 and /c«(2) - Ao(2) + e^/d = 10.9 
eV[l], yielding A = 2.92 eV. Besides setting the energy 
scale for gap excitations, A largely governs the magnetic 
exchange energy J. Using a canonical standard set of 
parameters in the limit of small hole hoppingj^, [l^], an 
exchange energy J ^ 300 meV is obtained. This value 
larger than that found in experiments (ll| results from 
the limitations of a perturbative expression for experi- 
mentally relevant parameter sets. 

The Ewald calculation was then repeated with a single 
oxygen dopant atom inserted into the unit cell, shown 
schematically in Fig. [1^, and the neighbouring atoms 
displaced as indicated from recent LDA studies [1^. The 
oxygen dopant was assigned formal valence, with surplus 
charge distributed equally among orbitals in the closest 
Cu02 planes. 

The site-dependent Madelung energies are presented 
in Fig. [l)D-d, showing large scale variations for sites 
closest to the dopant atom. Suppressions/enhancements 
of 0/Cu Madelung energies are observed, respectively, 
rising/falling to bulk values, shifted by the presence of 
doped holes, further away from the dopant. However, 



since the relative sign of the Madelung energies for Cu 
& O are negative, these changes largely cancel for A$a/ 
and thus A (Fig. [TJ;) is largely unaffected. 

Modifications in A allow us to examine the effect of 
the dopant on the exchange energy J (Fig. [Tf). Reflect- 
ing the spatial variations of A, J is suppressed by up to 
20 % near the dopant. Allocating dopant excess charge 
on either on Cu, on O, or distributed on both result in 
slightly different values of J (Fig. [T^) but the overall 
suppression seems rather immune to the way in which 
charge is distributed. 

To test this in a non-perturbative way ED studies of 
three-hole CU2O7 (Fig. 2a) and CuaOg (Figs. 2b & 2c) 
clusters were employed to determine changes to the ZRS 
Including Cu idr^2_ 
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orbitals, the Haniiltonian is H = ^ -ffi.o- where: 
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where i5 denotes the Cu-0 basis vectors and n the number 
operator. From this 3-band model, an effective single- 
band Hubbard or t — J model is derived, with effec- 
tive nearest 2t and next-nearest 2t' neighbour hoppings 
determined from the bonding/anti-bonding splitting of 
the ZRS, and J derived from singlet-triplet splitting. In 
terms of clusters CU2O7 determines t and J as the ZRS 
involves a common bridging oxygen while CU2O8 clusters 
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FIG. 2: (a) CU2O7 cluster used to extract t and J. (b) and 
(c) Two CU2O8 clusters shown in red outline used to extract 
t' and J'. The values for ep^d are indicated, (d) modifications 
to tpd due to the dopant induced displacements. 
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FIG. 3: (al), (bl) Normalized local crystal fields at the pla- 
nar and apical oxygen sites of the doped lattice. (a2), (b2) 
The corresponding electron phonon coupling strength \ cc . 
The data points follow the same key as figure 1. 



yield t' , J' via 0-0 hopping |13l|. For the undoped lattice, 
t = 383, t' = —141 meV and antiferromagnetic exchange 
couplings ofJ = 218, J' = 17 meV are characteristically 
obtained. J may be fine-tuned by adjusting tpd and tpp. 

The dopant is included by locally varying the on-site 
Cu and O energies ed and ep, respectively, shown in 
Fig. 121 Here two CU2O8 clusters are used which dif- 
fer with respect to the location of the dopant, and the 
site energies have been determined from site modified 
Madclung energies (Fig. 1), according to Eq. [TJ For 
the CU2O7 cluster t — 390, J = 216 mcV are obtained, 
while for the two CU2O8 clusters (Figs. |21d & |2t) we 
obtain t' = -242(-189) mcV and J' = 18(17), respec- 
tively. Electrostatic modifications to J and J' are slightly 
suppressed over undoped cluster values, in general agree- 
ment with Madelung estimates, although the magnitude 
is smaller. Importantly, we note that the symmetric 
placement of the dopant for the CU2O7 cluster gives only 
small changes to t, while for CU2O8 the increases are not- 
icably larger, particularly for geometry Fig. j^b compared 
to 2c. The asymmetric location of the dopant favors oc- 
cupation of the oxygen ligand orbitals in the plaquctte 
containing the dopant, giving larger modification of t' . 

We have repeated these cluster calculations, including 
the modulations in tpd induced by the structural distor- 
tions, according to the LDA displacements of Cu-0 bond- 
distances[l2. 14|. shown in Fig. 2d, and obtain t = 374, 
J = 197, t' = -231(-180) and J' = 17(16) meV. This 
tends to further suppress J. Thus oxygen dopant's net 
effect is to mildly suppress J and increase t' , indicat- 
ing that the dopant cannot be viewed as only modifying 
site energies and increasing J in downfolded single band 
models [7|. 

Since the Madelung energies are locally modified on 
the eV scale, real space modifications may result in sub- 



stantial changes to local crystal fields. Ions vibrating 
along the c— axis are sensitive to spatial gradients of the 
Madelung energies, and Ref . [l5| has shown that these val- 
ues of the local field determines the overall strength of 
cl-ph coupling at the oxygen sites, controlling coupling to 
Raman active Aig planar and apical oxygen vibrations, 
and out-of-phasc planar Big vibrations. Therefore, we 
have examined the dopant's effect on the crystal field in 
connection to local changes in el-ph coupling. 

For the undoped crystal, c-axis crystal fields are 
Eapex,piane = 18.74, 1.16 cV/A for the apical, planar oxy- 
gen sites, respectively, oriented as indicated in Fig. [T^. 
For the doped lattice, in Figs. [3^1, bl we plot the local 
c— axis electric field strength at the planar, apical oxy- 
gen sites, respectively, and the corresponding el-ph cou- 
pling strength A in Figs. [3^2, b2[l^. One can see that 
the i?-field strength is very sensitive to the local symme- 
try breaking introduced by the dopant, especially in the 
case of the planar oxygen atoms. The presence of the 
dopant's bare charge in the otherwise positively charged 
SrO/BiO structure suppresses the field in the closest ly- 
ing plane and the structural changes further modify the 
local fields. The largest changes to the i?-field occur in 
the plane whose field is oriented towards the dopant. In 
this case the geometry is such that the dopant's bare 
charge increases the strength of the original field, driving 
A for the c-axis planar oxygen modes up by a factor of 
5. As it is well known that the c— axis phonons give an 
attractive interaction in the d^^-y^ channel, this suggests 
that superconductivity may be locally promoted by the 
dopant, in agreement with the assessment of Ref. @. The 
enhanced el-ph coupling may on the otherhand drive a 
tendency to locally bind a hole to the lattice near the 
dopant location. This raises the possibility that LDOS 
modifications could be related to local polarons rather 
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FIG. 4: (a) The electron addition and removal spectra for 
the Cu02 cluster coupled to c-axis oxygen vibrations as a 
function of el-ph coupling strength. The energy gap Eg is 
indicated, (b) The electron addition (positive energy) and 
removal (negative) spectra as a function of energy for selected 
coupling strengths. (c,d) Eg, J as a function of el-ph coupling, 
respectively. 



than pairing. 

To address this issue, we examine three-site polarons 
via ED for an isolated Cu02 Hubbard cluster cou- 
pled to c-axis oxygen vibrations: H = H^i + Hiat + 
Hei~ph, where Hgi is defined in Eq. [2l and Hiat = 
^.yfi", Hei^ph = Y.u,s.a 9u{bl + )e^pJ_^Pd\<T • Here 
fi'^ is the phonon number operator for branch is the 
polarization of the ;^-th quantized local displacement, and 
cji, = eEy^h/2Mo^u sets the el-ph coupling strength to 
mode v. We consider coupling to uniform out-of-phase 
Big and in-phase Aig c-axis modes (fJ = 36 and 55 mcV, 
respectively) coupled to the oxygen hole density by the 
local field eE. We have not included any local modifi- 
cations to the phonon mode energies compared to the 
bulk. 

To make contact to STM LDOS, the electron addition, 
removal spectra A± , defined as 

A±(c.) = V |(v&^'|c, c^\^l)\'S{i. - E^'' + El), (3) 



where VE*" denotes the n-holc states with energy Ef, are 
plotted in Figs. 4a & 4b for different values of el-ph cou- 
pling gsig- A large number of phonon quanta have been 
used such that the results are independent of the num- 
ber of phonons for the parameter region investigated. As 
the coupling is increased, the spectral weight is gradu- 
ally transferred into phonon side bands giving broader 
and suppressed spectral peaks. The corresponding en- 
ergy gap Eg between the first removal and addition states 
decreases with increasing gsig (Fig. Ht) as the effective 
charge transfer energy is reduced by the gain in lattice en- 
ergy. As a consequence, the value of J (correspondingly 



determined from 2-hole CU2O7 clusters with phonons) 
increases substantially with increasing gsig, as shown in 
Fig. 4d. The decrease of Eg and increase of J due to 
local phonons will overwhelm the countering effects from 
purely electrostatic considerations as the coupling to the 
lattice increases. 

To estimate the size of the effect for Bi-2212, Epiane — 
1.16 eV for the undoped lattice yields gs^^ ^ 0.073 cV, 
well into the large polaron regime where side-bands are 
weak. Near the dopant however gsig is enhanced to ~ 0.2 
cV, where side bands begin to develop spectral weight in 
the removal/addition spectra and the charge gap is re- 
duced (Figs. [4^-c). At the same time, for this parameter 
regime, which is characteristic of multilayer cuprates, J 
is increased by 20 meV. much greater than the small re- 
duction determined from electrostatic effects alone. Thus 
for realistic parameter regimes, c-axis phonons act in con- 
cert with strong spatial variations of Madclung energies 
give an increase in J as well. The dopant may then pro- 
vide two coupled channels for d— wave pair enhancement, 
while causing suppressed spectral features as a conse- 
quence of strong local e-ph coupling. This is qualitatively 
what has been observed in experiments Q. Such synergy 
between phonons, polarons and antiferromagnetism has 
been already noticed in cluster quantum Monte Carlo 
studies [lit. 

In summary, we have investigated electrostatic modu- 
lations of local Madelung energies arising from the pres- 
ence of a dopant atom in Bi-2212 unit cells. While eV 
changes are found for the Madelung energies for copper 
and oxygen, these changes largely cancel for the charge 
transfer energy and give small local suppression to J. 
However the strong local variations in Madelung ener- 
gies are manifest in order 1 changes in el-ph coupling for 
c— axis oxygen modes. Using cluster studies, it was found 
that in combination electrostatic modifications and cou- 
pling to the lattice yield broadened spectral features, re- 
duced charge gap energies Eg , and a sizable local increase 
of J, implying a strong local interplay between anti- 
ferromagnetism, polarons, and superconducting pairing. 
The amount of variation of the local charge gap can thus 
serve as an important diagnostic for determining lattice 
coupling, electrostatic effects, and pairing. It is an open 
question whether a link between these quantities can be 
made experimentally. 
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